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Abstract. We present preliminary results from high resolution magneto-convection simulations 
where we find the formation of flux concentrations from an initially uniform magnetic field. We 
compute the effective magnetic pressure but find that the concentrations appear also in places 
where it is positive. The structures appear in roughly ten convective turnover times and live 
close to a turbulent diffusion time. The time scales are compatible with the negative effective 
magnetic pressure instability (NEMPI), although structure formation is not restricted to regions 
where the effective magnetic pressure is negative. 
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1. Introduction 

The current paradigm of sunspot formation is based on the idea of buoyant rise of 
flux tubes from the base of the solar convection zone to the surface of the Sun (ParkeiJ 
[l955„). This process is parameterise d in the widely used flux transport dynamo models 



(e.g. iDikpati fc Charbonneaij 1999[ ) in the form of a non-local a-effect where the strong 



(around 10^ G) toroidal magnetic fields in the tachocline give rise to poloidal fields at the 
surface. This poloidal magnetic field is then advected by meridional circulation back to 
the tachocline where it is amplified. 

These concepts face several theoretical difficulties, however, including the storage and 
gener ation of strong magnetic fields beneath the convection zone fe.g. Guerrero fc Kapvlgj 
2011 ). and the stability of the tachocline in the presence of such strong fields (jArlt et al.l 
20051). Furthermore, observations of sunspot rotation suggest that they might be a shallow 
phenomenon possibly occurring within the near surface shear layer (Brandenburg„.2005.,) . 
This requires a new mechanism to form sunspots. 

Theoretical works have sho wn that suitable turbulence can have a negative contribution 
to the magnetic pressure (e.g. iKleeorin et al. I I1990I Il996l : iKleeorin fc Rog achevskiilll994l 



Rogachevskii fc KleeorinI [20071 and references therein). This effect leads to the negative 



effective magnetic pressure instability (NEMPI) where even uniform, sub-equipartition, 
magnetic fields can form flux concentrations. This is compatible in view of the results 
from direct simulations (DNS) of convection driven dynamos where diffuse magnetic 
fields are generated in all of the convection zone (e.g. iGhizaru et al. 2010l : Kapvla et al 
2012al ). 



Recently, a lot of effort has been devoted to study this effe c t ^sii i g mean-field model s 
and DNS of forced turbulence (e.g. ' Brandenburg et al."2010' "2012; Kemel et al. "20121), 
culminating in the detection of NEMPI in DNS (Brandenburg et al. 2011). A negative 
turbulent contribution to the effective (mean-field) magnetic pressure has also been found 
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for convection (^Kapvla et al. l l2Q12bh but no NEMPI has been detected so far. Here we 



present results from new high resolution convection simulations designed to be better 
suited for the detection of NEMPI. 



2. Model 

We solve the compressible hydromagnetics equations, 

dA 

— = ux B - r]/j.oJ, (2.1) 

^ = ^ + - (V • 2iypS -Vp^JxB), (2.3) 
JJt p 

T— = - [V • {KVT + XiPTVs) + ii^r^J^] + 2vS', (2.4) 
Vt p 

where A is the magnetic vector potential, u is the velocity, S = V x A is the magnetic 
field, J = Pq^^ X B is the current density, r] is the magnetic diffusivity, po is the vacuum 
permeability, D/Dt = d/dt -\- u • ^ is the advective time derivative, g = —gCz = const 
is gravity, v is the kinematic viscosity, K is the radiative heat conductivity, Xt is the 
unresolved turbulent heat conductivity, p is the density, s is the specific entropy, T is the 
temperature, and p is the pressure. The fluid obeys the ideal gas law with p = (7 — l)pe, 
where 7 = j c\j = 5/3 is the ratio of specific heats at constant pressure and volume, 
respectively, and e = cyT is the internal energy. The traceless rate of strain tensor S is 
given by 

We omit stably stratified layers above and below the convection zone. The depth of 
the layer is Lz = d whereas the horizontal extents are Lx/d = 10 and Ly/d = 5. 
The boundary conditions for the flow are impenetrable and stress free, and perfectly 
conducting for the magnetic field. The energy flux at the lower boundary is fixed, and 
we use a black body boundary condition given by 

aT^ = -CPPX^ - pTXt^, (2.6) 



where a is the Stefan-Boltzmann constant at the surface (cf. iKapvla et al. l l2Qllh . 



We use a constant x = K/{cy>p) whereas Xt is zero below z/d < 0.5, Xt/x = 10 i^i the 
range 0.5 < z/d < 0.9, and Xt/x = 50 above z/d> 0.9. The Prandtl number Vi = v/x is 
equal to 10. In this setup convection transports the majority of the flux whereas radiative 
diffusion is only important near the bottom of the domain. We start a hydrodynamic 
progenitor run from an isentropic stratification with density stratification of 80. The 
density and pressure scale heights, mean entropy profile, equipartition magnetic field 
5eq = {PQpu^)^^'^ ^ and the Mach number. Ma = iirms/cs, in the thermally saturated 
state of the simulation are shown in Fig. [TJ In the simulation considered here the fluid 
and magnetic Reynolds numbers are Re = 'Urms/(^^i) ^ 94 and Rm = Urms/iv^i) ^ 5, 
respectively, where Urms is the rms value of the volume averaged velocity and ki =27r/d. 
We use a grid resolution of 1024 x 512 x 256. The computations were performed with the 
Pencil CodeQ. 



t http:/ /pencil-code. googlecode.com/ 
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Figure 1. Left panel: specific entropy (solid line), and density (dashed) and pressure (dot- 
dashed) scale heights as functions of depth. Right panel: Mach number (solid line) and equipar- 
tition field strength from the thermally saturated regime as functions of depth. 




Figure 2. Vertical velocity Uz near the surface (left panel) and near the bottom (right) of the 
domain from a hydrodynamical simulation. The contour levels correspond to zb0.37Ma (left) 
and zb0.03Ma (right). 



3. Results 

We first allow a hydrodynamic progenitor simulation to saturate after which we im- 
pose a uniform horizontal field BoCy with Bq/Bqc^ ^ 0.3. The vertical velocity Uz near 
the surface {z/d = 0.98) and near the bottom (z/d = 0.02) of the domain in the hy- 
drodynamic run are shown in Fig. [2l In forced turbulence simulations NEMPI appears 
wh en the scale separation b etween the forcing scale and the box size is of the order of 
15 ([Brandenburg et al.ll201l[ ). In our convection setup this is probably satisfied near the 
surface but not near the bottom of the domain. 

We find that large-scale structures form within ten convective turnover times; see Fig. [3] 
where By = By{x^z) — {By)xy{z) is shown. The subscripts refer to averages over either 
only y or over both horizontal directions, respectively. We remove the horizontal mean 
value (see the left panel of FigH]) in order to make the horizontal variation of Byix^z) 
visible. The magnetic structures appear near the surface and sink on a timescale of a 
few tens of turnover times r = (i^rms^f where we estimate kf by = cjrms/'^rms, 
and (jj = V X n. Now the turbulent diffusion time Tdiff = (^t^i)~^ is roughly 180 
turnover times if we assume r]t = Urms/ (3A:f ) for the turbulent diffusivity. The maximum 
field strength in the concentrations is of the order of the imposed field. The elongated 
magnetic structures are also weakly discernible in the instantaneous magnetic field By^ 
see the left panel of Fig. [5l In the vertical field B^ from the same instant, however, it is 
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not possible to distinguish the same structures, right panel of Fig. O This can due to the 
perfect conductor boundary condition which imposes 5^ = a the boundary. 

We define the effective magnetic pressure as T^eff = ^ where Qp = —2/S.pui — 

eq " 

A6^ + 2Aby^ and where A refers to the difference between runs with and without an 
imposed field. In the present case a small-scale dynamo is absent and thus the magnetic 
correlations come only from the run with the imposed field. We find that the effective 
magnetic pressure T^eff is positive near the upper boundary and negative below z/d < 0.1 
with increasingly negative values towards the bottom, see the right panel of Fig. IH The 
maxima of T^eff are associated with maxima of By, whereas the minima of T^eff coincide 
with lower average values of By, see Fig. [3l 



4. Conclusions 

Our results have shown that convection can lead to magnetic fiux concentrations by 
a mechanism that may be related to NEMPI. A possible alternative mechanism is mag- 
netic fiux expulsion that has previously been found to be responsible for a segregation of 
magnetised and unmagnetised regions in large aspect ratio convection simulations with 
an imposed vertical magnetic field (T ao et al. 1998) . Our results appear to be similar, 
except that here we have an imposed horizontal mag:neti c field. In that respect, it is 
useful to mention recent simulations of IStein fc Nord lund' (2012), who inject a 1000 G 
horizontal magnetic field at the bottom of their simulated convection domain and find 
after some time the emergence of a bipolar magnetic field at the surface. In their case 
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Figure 4. Left panel: horizontally averaged horizontal magnetic fields Bx and By in units of 
the imposed field Bq. Right panel: horizontal average of the effective magnetic pressure T^eff- 
The average is taken over a period where the rms magnetic field is saturated. 




Figure 5. Magnetic field components By (left panel) and Bz (right) near the surface 
{z/d = 0.98) from the same time as in the bottom panel of Fig. O The contour levels cor- 
respond to -Bo < By - Bo < Bo (left) and -Bo < B^ < Bo (right). 



the reason for the formation of flux concentrations is argued to be the downdrafts of 
the deeper supergranulation pattern, which tend to keep the magnetic field concentrated 
into flux bundles at the bottom of their open domain. Our present simulations do not 
capture this effect, because they are probably not deep enough and our domain is impen- 
etrative and perfectly conducting at the bottom, excluding therefore their mechanism as 
a possible explanation. Of course, another important difference between our simulations 
and those of ilStein fc Nordlundl (|2Q12l) is the presence of a radiating surface in their case. 
This might enhance magnetic flux concentration s formed throujSi:h loc al suppression of 
convective energy flux by magnetic fields (Kitchat inov fc Mazurll2QQQ[ ). This might well 
be an important effect that needs to be studied more thoroughly. 
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Figure 6. Effective magnetic pressure T^eff from the same times as in Fig. [S] 
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